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Abstract: Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. 
The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine 
time resolution spectral analysis shows power-law decays of the peak energy from the onset of 


the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly 
correlated power-law behavior is observed between the luminosity and the spectral peak energy 
that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any 
of the existing models to account for all of the observed spectral and temporal behaviors 
simultaneously. 

Main Text: In the context of gamma-ray bursts, GRB 130427A, which triggered the Gamma-ray 
Burst Monitor (GBM) (7) on the Fermi Gamma-ray Space Telescope on 2013-Apr-27 at To = 
07:47:06.42 UTC (2 - 4) is an extreme case. The peak flux on the 64 ms timescale is 1300 ±100 
photons s 1 cm 2 in the 10 - 1000 keV range and the fluence, integrated over the same energy 
range and a total duration of approximately 350 s, is (2.4 ± 0.1) x 10" 3 erg cm 2 . The longest 
continuously running GRB detector, Konus on the Wind spacecraft, has been observing the 
entire sky for nearly 18 years and only one burst had a larger peak flux, by -30% (GRB 
1 109 18 A) (5). GRB 130427A is the most fluent burst in the era starting with the 1991 launch of 
the Burst And Transient Source Experiment (BATSE) on the Compton Gamma-Ray 
Observatory. Finally, the energy of the spectral peak in the first time bin (To -0. 1 to 0.0 s), 5400 
± 1500 keV, is the second highest ever recorded (6). 

The initial pulse (Fig. 1), lasting up to 2.5 s after the trigger, stands on its own as being so 
bright (170 ± 10 ph s -1 cm 2 peak flux for 10 - 1000 keV in the 64 ms time bin at To +0.51 s) as 
to be ranked among the 10 brightest GBM or BATSE bursts (7-9). The brightness allows us to 
track the spectral evolution of the rising portion of a well-separated pulse with unprecedented 
detail (10). Evident in the GBM low-energy light curve (Fig. 1; as well as the 15 - 350 keV light 
curve presented in (11)) are fluctuations starting at around 1 s that are not present at higher 
energies. If these represent additional low-energy pulses, their presence clearly does not 
dominate the analyses presented below. 

Past studies of time-resolved spectra of simple pulses in GRBs indicate that there are 
broadly two classes of spectral evolution. These are called ‘hard-to-soft’ and ‘tracking’ pulses 
(12, 13), depending on whether the energy of the peak in the vF v spectrum (generically called 
Apeak herein) monotonically decays independently of the flux evolution or else generally follows 
the rise and fall of the flux. Typically, there are at most one or two spectra available for fitting 
during the rising portion of the flux history. What makes this event unique is that there are 
roughly 6 time bins with excellent counts statistics before the peak in the 10 - 1000 keV flux. 

As seen in Fig. 1, there is a clear trend in the individual detector’s light curves: the > 20 
MeV Fenni Large Area Telescope (LAT) low-energy (LLE) (14, 15) light curve peaks before the 
GBM trigger time (Jo), while the GBM bismuth germanate (BGO) detector #1 (300 keV - 45 
MeV) and sodium iodide (Nal) detector #6 (8 - 300 keV) peak at successively later times. To 
quantify this, we performed an energy-dependent pulse lag analysis using a Discrete Cross 
Correlation Function (DCCF) and obtain the time lags r (16) between the highest energy LAT 
LLE light curve and light curves at several selected energy ranges in the GBM Nal and BGO 
detectors (Fig. 1 - inset). We find good agreement between the expected lag behavior and the 
pulse width model W(E) oc E a (17), obtaining a fitted value for a = -0.27 ± 0.03 (18). This 
model was previously fit to 400 pulses from 41 BATSE GRBs (17); an average value of a = - 
0.41 was found. Synchrotron shock model simulations made by (19) found a — 0.4 for pulses 
of 2 - 10 s duration but a > — 0.2 for pulses of 0. 1 - 1 s. Three LAT photons with energies 



greater than 100 MeV are clustered in coincidence with the LLE peak, and so may arise by the 
same mechanism. 

Although most GRB spectra are well fit by the smoothly joined broken power-law 
function of Band et al. (20), in some cases the simultaneous fit of a Band function together with 
an additional blackbody component is significantly better statistically (6, 21, 22). Burgess et al. 
(23, 24) also noted the requirement of an additional blackbody component, but replaced the 
phenomenological Band function with a physically-motivated synchrotron function. We present 
two separate time-resolved spectral analyses for the first 2.5 s of GRB 130427A, with 
comparable goodness of fit: the Burgess et al. synchrotron function plus blackbody, and the Band 
function (18). A blackbody component is not required when the more flexible Band function 
alone is used. Although the time evolution of £ pea k as determined by Band function fits is 
consistent with a single power-law (with an index of -0.96 ± 0.02), the evolution of the 
synchrotron peak energy is not (Fig. 2). A broken power-law fit is better constrained and shows a 
shallower decay before the pulse peak, with an index of -0.4 ± 0.2 during the rising phase and - 
1.17 ± 0.05 during the decaying phase and a fitted break time at T 0 + 0.28 ± 0.08 s, or ~ 0.2 s 
before the pulse peak in the 10 - 1000 keV flux. Both fitted indices during the decay phase are 
consistent with the -1 power law index expected from standard fireball curvature effects (25, 

26). The shallower spectral peak decay index prior to the light curve decay phase has a natural 
explanation in the context of the pulse being driven by a shock between thick colliding shells 
(19). 

Two-component models including a thennal contribution (27) constrain the value of the 
photospheric radius using the blackbody flux and temperature (kT - see Table SI). Comparing 
with the flux of the dominant non-thermal spectral component then permits determination of the 
Lorentz factor at the photosphere (T p h) (28). As shown in Fig. 3, the minimum value of the 
photospheric bulk Forentz factor /j,h starts out at 500 and monotonically decreases to ~ 1 00 over 
the duration of the pulse (similar to behavior observed in GRB 1 10721 A) (29). Internal shocks 
require higher Forentz factors at later times. However, this might still be consistent with the 
monotonically decreasing 7j,h if the outflowing shell that produces this photospheric component 
produced the non-thermal triggering pulse by colliding with a slower and slightly earlier ejected 
shell that did not produce detectable photospheric emission. Otherwise, the observed behavior 
would favor magnetic reconnection models or mini -jets (30, 31), which abandon a simple 
spherical geometry. 

Using the measured redshift of z = 0.34 (32), the host rest-frame luminosity and 
synchrotron peak energies are calculated and the decay phase apparent isotropic luminosity L- 
Fpeak correlation is fit with a power-law index of 1.43 ± 0.04 (Fig. 4). A theoretical analysis of 
high-latitude curvature radiation produced in relativistic shell collisions of spherical blast waves 
shows that L oc Fj, eak 3 during the decay phase of a pulse (25, 26), contrary to the behavior shown 
in Fig. 4. In a picture of an expanding fluid element rather than a colliding shell, synchrotron 
emission by electrons with characteristic energy y e obeys the relations E pea k °c / "By 2 and L oc 

T 2 B 2 y 2 , for 7c and B both in the jet frame. In the optically-thin coasting phase of the outflow, the 
bulk Forentz factor r is constant. Naively assuming that the magnetic flux is frozen in the flow 
(BR 2 oc const, in the co-moving frame, where R is the comoving emission region radius), then 
adiabatic losses of the electrons imply y e oc R A short calculation then gives L oc E pe a k 3/2 , 
which is consistent with the 1 .43 index derived from the data. A constant expansion velocity 



dR/dt scenario predicts, however, f pca k R 4 oc t 4 . Other jet-wind assumptions yield different 
correlations, e.g., for the deceleration epoch where dT/dt < 0, or for radial field evolution 
appropriate for jet cores. 

The isolated initial pulse of GRB 130427A is apparently unmodified by preceding engine 
activity or nascent external shock emission. Our analysis shows that there is good agreement 
between the pulse width as a function of energy and the expected lag, the characteristic energy 
has roughly a -1 power-law decay with time during the decaying phase, the temperature of the 
blackbody component implies a photospheric radius that is incompatible with the internal shock 
radius, and the apparent isotropic luminosity is related to the 3/2 power of the characteristic 
energy. It is a challenge to explain all these behaviors simultaneously. 
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Fig. 1 . The first 3 s of GRB 130427A. Shown are composite light curves for the three 
Fermi detector types {green: GBM Nal #6 [10 - 300 keV]; blue: GBM BGO #1 [300 
keV - 45 MeV]; red: LAT LLE [> 20 MeV]). Each curve has been normalized so that 
their peak intensities match. High probability LAT photons >100 MeV are indicated by 
circles {right axis - energy in MeV). {Inset Figure:) Lag analysis of the triggering pulse 
of GRB 130427A. Time lag r {filled symbols) as determined by the DCCF analysis 
between the (10-100 MeV) LLE lightcurve and selected energy bands of the Nal {green) 
and BGO {blue) lightcurves. Also displayed are fitted pulse widths as a function of 
energy W(E) ( hollow symbols, in sec.) for several energy bands. The two dashed lines 
represent: 1), the best-fit power-law model (y 2 of 5.6 for 9 degrees of freedom) for W{E) 
{black), and 2), the expected dependence of the time lag rasa function of energy {red), 
assuming the same power-law index as in 1). 



Time [s] 

Fig. 2. The fitted Band function A pea k {blue) and synchrotron peak energies (red) as a function of 
time. The times are referenced from when the LLE light curve peaks 0.1 s before the trigger. A 
broken power-law fit to the red points is indicated by a dashed line (early time decay index of- 
0.4 ± 0.2, with a break at 0.38 ± 0.08 s, breaking to an index of-1.17 ± 0.05 with a "i = 28 for 22 
degrees of freedom). We also show the Band function Apeak values for the same time intervals 
with a single fitted power law index (-0.96 ± 0.02 with a y 2 of 1 9 for 24 degrees of freedom). 
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Fig. 3. Plot showing trends in the derived photospheric Lorentz factor ( red - left axis ) and radius 
{blue - right axis). The reference time is the same as in Fig. 2. We obtain both values from the 
instantaneous ratio of the observed blackbody component flux to the total flux, following Eq. 4 
& 5 in (25) and assuming a value of To = 1 for the ratio between the total emitted thennal energy 
vs. the total energy emitted in gamma rays. 


[cm] 



Rest Frame E peak [keV] 

Fig. 4. Correlation between the GRB 13 042 7 A host rest frame synchrotron peak energy and 
isotropic luminosity during the rising phase ( red) and decaying phase {black) of the triggering 
pulse. Time progresses approximately from right to left on the plot. The 1.43 ± 0.04 power-law 
index fit to the black points is shown in green (region of uncertainty in grey), while the 3/2 
power-law from the magnetic flux-freezing calculation in the text is indicated in blue. 
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Spectral Analysis Method: Band Function Example 

The GBM time-tagged event (TTE) data (1) and the LAT low-energy events (LLE) data (14, 15) 
were binned at 100 ms resolution, referenced from the GBM trigger, covering a continuous 
energy range from 8 keV to 100 MeV. We used GBM data for BGO #1, plus Nal #6, with 
source-to-zenith angles of 53.5° and 6.5° respectively, as determined for the best available 
CARMA radio source position (173.1367° RA, +27.6989° Dec ± 0.4", J2000) (33). Using the 
spectral analysis package RMFIT (34), we first fit the Band GRB function (20) to the entire 
interval from To -0. 1 to 2.5 s in order to determine the relative offset between the BGO and Nal 
rates, since there is typically an effective area correction of 10% between the BGO and any of 
the Nal detectors. We fix the effective area correction between the BGO and LLE rates to 1.0, as 
there is no known offset between the two and unfreezing the parameter drives it to large, 
unreasonable values. The best-fit Band parameter values are fs pea k = 376 ± 7 keV, as (the low- 
energy power law index) = -0.86 ± 0.01 and fi (the high-energy power law index) = -2.66 ± 
0.035, with C-Stat = 599.65 for 377 degrees of freedom and poor residuals. With this spectrum 
as a seed, we then perform spectral fits for each of the time bins in the series. Examination of the 
results indicates that the /? parameter is undetennined in roughly half the spectra and when it 
could be determined by the fits, /? never deviated beyond one standard deviation from the value 
detennined from the fit to the entire interval. We find ) 3 to be consistent with a constant value of 
-2.66 and so the results shown are based on fixing /? to that value and fitting the other 
parameters. The spectral evolution of the £p ea k parameter is evident in Fig. 2, falling smoothly 
from the onset of the pulse; the complete set of spectral fit parameters can be found in Table SI. 
Finally, the distribution in an peaks very close to -0.66, or -2/3, the value expected for optically 
thin synchrotron emission. 

If one allows for the evolution of a magnetic field in a fast-cooling synchrotron scenario, low- 
energy spectral indices of -0.8 are possible, as shown by (35). The fitted an values during the 
brightest portion of the pulse fall in the -0.2 to -0.6 range, many sigma away from the desired - 
0.8 spectral index. It is usually the case that if one then imposes a low energy spectral index 
value (say, -0.8), a second component, such as a blackbody, is required to make up the 
difference. So, we would have to allow a blackbody component along with the quasi-Band 
spectrum, in this purely magnetic model. Interestingly, the Band an does settle on an average that 
is nearly -0.8 during the later decay phase. 

Lag Analysis 

An empirically-motivated function for the expected time lag between two pulses of energy E\ 
and E 2 is given by t(E\, E 2 ) = Use (E\ a - E 2 a ), where their full widths at half maximum (FWHMs 
- in s) are related by the function W(E) = Wo (E/E 0 ) a , where E 0 = 1 keV (17). We fitted pulses at 
several energies with a lognormal function to obtain the FWHM in each energy band (Fig. 1 - 
inset, hollow symbols). These were then used to derive the parameters Wo = 3.2 ± 0.5 s and a = - 
0.27 ± 0.03. We find good agreement between the above expected lag behavior, using our fitted 
value for a (Fig. 1 - inset, dashed lines) and the values of t V]SC determined by the width fitting. 

Photospheric Radius and Lorentz Factor 

Given the joint observations of the photospheric radius and Lorentz factor, how consistent are 
these with the standard internal shock model of GRBs? The radius ris of the corresponding 



internal shell collisions then should be ~ c T 2 At ~ 7.5 x 10 15 (77500) 2 cm » r p h~ 10 12 cm, 
where At ~ 1 s is the pulse duration (see Fig. 1), which clearly is not the case, as the 
photospheric and synchrotron spectral components are so closely coupled in time. Such large 
radii are also required for the y-ray emission region to be optically thin to yy absorption, which 
would cut off the spectrum at lower energies than observed. Because the wind escaping the 
photosphere travels more slowly than light, the emission from the colliding shells would lag 
behind the photospheric radiation by -At. Finally, we note that r p i, and /"scale with the ratio 
between the total gamma-ray energy and the thennal energy Yo, as described in Fig. 4. Thus, r ~ 
500 should be considered a lower limit. A more physical value for To must be greater than 1, 
making r correspondingly larger and thus increasing the total isotropic energy of the burst at all 
energies. As shown in the Fermi LAT companion paper (10), the best estimate for the total 
apparent gamma-ray energy is ~10 34 erg, which is uncomfortably large if the total bolometric 
energy inferred from r is more than a factor o 1' To ~ 1 0 larger. 
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Time 

Range 

(s) 

Band £ pcak 
(keV) 

Band a B 

Band 

C-Stat 

(371 

dof) 

Synchrotron 

T peak 

(keV) 

Electron 
Index <5 

Photospheric 
Energy kT 

(keV) 

Total Flux 10 Delta C-Stat 

keV to 100 0 j 

Synch vs. 

Synch ±BB 

(y cnT 2 s~') 

Synchrotron 

C-Stat 

(/dof) 

-0.1 : 

: 0.0 

5400 ± 1500 

-0.9 ±0.1 

330 

2560 ± 1210 

3.8 ±0.5 

(no constraint) 

9 ± 1 

-0.8 

348 

/ 383 

0.0 : 

0.1 

1470 ± 170 

-0.49 ± 0.07 

389 

1800 ±590 

4.6 ±0.6 

276 ±34 

31 ±3 

39.1 

433 

/ 381 

0.1 : 

0.2 

849 ±58 

-0.25 ± 0.06 

365 

1420 ±320 

4.8 ±0.5 

175 ± 10 

46 ±4 

164. 

410 

/ 381 

0.2 : 

0.3 

687 ±38 

-0.33 ± 0.05 

358 

1270 ±220 

6.7 ± 1.5 

146 ±9 

85 ±5 

130. 

398 

/ 381 

0.3 : 

0.4 

553 ±29 

-0.46 ± 0.04 

369 

975 ± 98 

10.00 

129 ±9 

116 ± 6 

90. 

374 

/ 382 

0.4 : 

0.5 

504 ± 26 

-0.57 ±0.04 

396 

856 ±60 

10.00 

86 ±9 

154 ±6 

43. 

395 

/ 382 

0.5 : 

0.6 

380 ± 19 

-0.57 ±0.04 

364 

640 ± 44 

10.00 

70 ±6 

152 ±6 

44. 

361 

/ 382 

0.6 : 

0.7 

359 ±20 

-0.64 ± 0.05 

362 

587 ±39 

10.00 

54 ±9 

149 ±6 

12. 

376 

/ 382 

0.7 : 

0.8 

327 ± 19 

-0.73 ± 0.05 

362 

560 ±41 

10.00 

34 ±5 

137 ± 7 

19. 

411 

/ 382 

0.8 : 

0.9 

280 ± 18 

-0.76 ± 0.05 

396 

386 ±30 

10.00 

53 ±9 

126 ±6 

22. 

390 

/ 382 

0.9 : 

1.0 

253 ± 17 

-0.72 ± 0.06 

370 

338 ±30 

10.00 

49 ± 11 

106 ±6 

7.0 

365 

/ 382 

1.0 : 

1.1 

208 ± 16 

-0.76 ± 0.07 

357 

304 ± 25 

10.00 

30 ± 15 

95 ±6 

2.0 

380 

/ 382 

1.1 : 

1.2 

207 ± 16 

-0.69 ± 0.08 

327 

347 ± 40 

10.00 

22 ±3 

74 ±7 

1.2 

356 

/ 382 

1.2 : 

1.3 

182 ± 13 

-0.73 ± 0.08 

315 

238 ±27 

10.00 

43 ±8 

84 ±5 

20. 

350 

/ 382 

1.3 : 

1.4 

170 ± 11 

-0.62 ± 0.09 

393 

273 ±26 

10.00 

29 ±4 

85 ±6 

5.9 

396 

/ 382 

1.4 : 

1.5 

167 ± 12 

-0.69 ± 0.09 

363 

254 ±26 

10.00 

24 ±4 

79 ±7 

12. 

380 

/ 382 

1.5 : 

1.6 

165 ± 13 

-0.81 ±0.09 

342 

233 ±24 

10.00 

28 ±5 

73 ±6 

15. 

354 

/ 382 

1.6 : 

1.7 

130 ± 9 

-0.5 ±0.1 

363 

169 ±20 

10.00 

31 ±6 

73 ±5 

13. 

356 

/ 382 

1.7 : 

1.8 

142 ±11 

-0.7 ±0.1 

349 

193 ±20 

10.00 

29 ±5 

77 ±6 

7.2 

380 

/ 382 

1.8 : 

1.9 

128 ± 10 

-0.5 ±0.1 

297 

184 ±20 

10.00 

27 ±5 

69 ±6 

7.6 

321 

/ 382 



1.9 : 2.0 

137 ± 11 

- 0.7 ± 0.1 

339 

199 ±22 

10.00 

23 ±4 

68 ±7 

6.7 

346 / 382 

2.0 : 2.1 

120 ± 10 

- 0.7 ± 0.1 

341 

198 ±23 

10.00 

21 ±4 

60 ±7 

- 0.8 

317/382 

2.1 : 2.2 

114 ± 10 

- 0.6 ± 0.1 

366 

136 ± 21 

10.00 

27 ±4 

51 ±5 

18 . 

364 / 382 

2.2 : 2.3 

109 ± 11 

- 0.7 ± 0.15 

305 

221 ±41 

10.00 

17 ± 2 

41 ±6 

5.4 

300 / 382 

2.3 : 2.4 

109 ± 11 

- 0.8 ± 0.2 

299 

1 12 ± 32 

10.00 

35 ±8 

46 ±5 

8.1 

321/382 

2.4 : 2.5 

101 ± 11 

- 0.7 ± 0.2 

330 

170 ±39 

10.00 

13 ±2 

35 ±7 

- 2.4 

358 / 382 


Table SI. Spectral fit results for the onset pulse of GRB 130427A. Columns 2-4 contain the varying Band fit parameters ([:> was fixed) 
and C-Stat fit statistic, while the rest of the columns are derived from the synchrotron plus blackbody function spectral fit. The 
synchrotron peak energy (column 5) has not been scaled by the square of the electron minimum Lorentz factor. After 0.3 s, the 
electron distribution power law index (column 6) could not be constrained by the data and was thus fixed to the steep value of 10 
(essentially equivalent to thermal electrons). Column 9 indicates the change in the fitting merit function between the synchrotron 
model alone versus the synchrotron model plus blackbody. Although simulations are required to determine probabilities accurately, a 
value greater than 20 indicates that the extra blackbody component is required at the 4 sigma level, assuming normal statistics. 
Parameter uncertainties are 1 sigma statistical only. Times are relative to trigger time To. 




